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Abstract 

Shot noise affects differently the nonhnear electron transport in semiconductor 
superlattices depending on the strength of the coupling among the superlattice 
quantum wells. Strongly coupled superlattices can be described by a miniband 
Boltzmann-Langevin equation from which a stochastic drift-diffusion equation is 
derived by means of a consistent Chapman-Enskog method. Similarly, shot noise in 
weakly coupled, highly doped semiconductor superlattices is described by a stochas- 
tic discrete drift-diffusion model. The current-voltage characteristics of the corre- 
sponding deterministic model consist of a number of stable branches corresponding 
to electric field profiles displaying two domains separated by a domain wall. If the 
initial state corresponds to a voltage on the middle of a stable branch and is sud- 
denly switched to a final voltage corresponding to the next branch, the domains 
relocate after a certain delay time, called relocation time. The possible scalings of 
this mean relocation time are discussed using bifurcation theory and the classical 
results for escape of a Brownian particle from a potential well. 



1 Introduction 



Nonlinear cliarge transport in semiconductor superlattices has been the object 
of many theoretical and experimental studies in the past decade [1,2]. Super- 
lattices are artificial spatially periodic structures first proposed by Esaki and 
Tsu in order to realize a device that exhibits Bloch oscillations [3]. A su- 
perlattice (SL) in its simplest form contains a large number of periods, each 
comprising two layers, which are semiconductors or insulators with different 
energy gaps, but with similar lattice constants e. g., GaAs and AlAs. These 
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SLs are synthesized by molecular-beam epitaxy or related epitaxial growth 
techniques in the vertical direction. The conduction band edge of an infinitely 
long ideal SL is modulated so that in the vertical direction it looks like a 
one-dimensional (ID) crystal, which is formed by a periodic succession of a 
quantum well (GaAs) and a barrier (AlAs). Typical experiments of vertical 
charge transport use an undoped or doped SL of finite length placed in the 
central part of a diode (forming a p-i-n or n'^-n-n'^ structure) with respective 
contacts at either end of the diode. Depending on the bias condition, the SL 
configuration, the doping density, the temperature or other control parame- 
ters, the current through the SL and the electric field distribution inside the 
SL display a great variety of nonlinear phenomena such as pattern formation, 
current self-oscillations, and chaotic behavior [1,2]. 

To describe and understand nonlinear charge transport in SLs, it is essential to 
distinguish between weakly coupled and strongly coupled SLs. Weakly coupled 
SLs (WCSLs) contain rather thick barriers separating the SL quantum wells, 
i. e., the barrier width is much larger than the typical electron wavelength in- 
side the barrier. Therefore, a description of the electronic properties of WCSLs 
can be based on the subband structure of the corresponding isolated quantum 
well together with resonant tunneling across the barrier of two adjacent wells. 
In contrast, the quantum wells of strongly coupled SLs (SCSLs) are separated 
by thin barriers so that the electronic properties of SCSLs can be described in 
terms of extended states such as Bloch functions. The simplest mathematical 
models applied to a SL give rise to balance equations involving mesoscopic 
quantities such as the electric field, the electron density, the drift velocity, 
etc. The task of deriving these equations from first principles is far from be- 
ing completed, although reasonable deterministic balance equations have been 
derived for both SCSLs [4] and for WCSLs [2] in particular limiting cases. 

A fundamental difference between WCSLs and SCSLs is that the former are 
governed by spatially discrete balance equations (differential-difference equa- 
tions), whereas the latter are governed by spatially continuous equations (par- 
tial differential equations). Both types of equations may have solutions, whose 
electric field profiles display regions of high electric field coexisting with re- 
gions of low electric field. The resulting dynamical behavior is very different 
for these two types of equations. For SCSLs, which are described by continuous 
balance equations, the field profile consists of a charge dipole moving with the 
fiow of electrons, which very much resembles the Gunn effect in bulk semicon- 
ductors [5]. Under dc voltage bias, this basic motion results in self-sustained 
oscillations of the current through the SL due to the periodic movement of 
dipole domains. In contrast, in WCSLs, which are described by discrete bal- 
ance equations [6,2], electric-field domains (EFDs) are separated by a domain 
wall, which consists of a charge monopole. The domain wall in a WCSL may 
move with or opposite to the electron fiow, or is pinned, depending on the 
value of the current [7]. This pinning of the domain wall occurs only in the 
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discrete models. 



In this paper, wc are interested in modeling the effects of shot noise in SLs, 
and then discussing its influence on the motion of EFDs. Shot noise is a con- 
sequence of the quantization of the charge [8] . For SCSLs, we shall describe in 
Section 2 shot-noise effects by using a miniband Boltzmann-Langevin equa- 
tion, from which we shall derive a stochastic drift- diffusion equation in the 
hyperbolic limit by means of a consistent Chapman-Enskog method. The re- 
sulting stochastic drift-diffusion equation is similar to that describing fluctu- 
ations in the Gunn effect, and therefore enhancement of charge fluctuations 
near the threshold to self-sustained oscillations of the current is to be expected 
[9]. 



The discrete balance equations describing electron transport in WCSLs have 
not been derived consistently from a kinetic theory, and therefore we cannot 
follow a perturbative treatment as in the case of the SCSLs. Shot-noise effects 
will therefore be described in Section 3 by Langevin equations resulting from 
adding appropriate white-noise sources to the discrete equations [10]. An im- 
portant situation to observe the effects of shot noise is EFD relocation due to 
voltage switching [11]. The current- voltage characteristics of a WCSL consists 
of a number of stable branches corresponding to electric field profiles display- 
ing two domains separated by a domain wall. If the initial state corresponds to 
a voltage on the middle of a stable branch and is suddenly switched to a final 
voltage corresponding to the next branch, the domains relocate after a certain 
delay time, called relocation time. The deterministic theory of domain relo- 
cation is based on numerical simulations of a nonlinear differential-difference 
model [12]. One of the main effects of shot noise is to render the relocation 
time random. If the final voltage after switching is near the limit point that 
marks the end of the initial stable solution branch, the mean relocation time 
can be investigated using bifurcation theory. We calculate a possible scaling of 
this mean relocation time using a stochastic amplitude equation correspond- 
ing to the normal form of a saddle-node plus a term due to projected shot 
noise. For this equation, the classical results for escape of a Brownian particle 
from a potential well yield the scaling of the mean relocation time. The last 
Section contains our conclusions. 



2 Shot-noise effects in strongly-coupled superlattices 



A deterministic simple model of one- dimensional (ID) electron transport in a 
SCSL consists of the following Boltzmann-Poisson system [4] : 
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Here Z, £, /, n, A^c, /c^, T, — F, m* and — e < are the SL period, the dielec- 
tric constant, the onc-particle distribution function, the 2D electron density, 
the 2D doping density, the Boltzmann constant, the lattice temperature, the 
electric field, the effective mass of the electron, and the electron charge, re- 
spectively. The right side of Eq. (1) contains two collision terms. The first one 
has BGK (Bhatnagar-Gross-Krook) form, and it represents inelastic energy 
relaxation towards a ID effective Fermi-Dirac distribution /™(/c;n) (local 
equilibrium) with collision frequency z/g. In Eq. (4), the chemical potential fi 
depends on n and is found by inverting the exact relation (3). The second 
collision term accounts for impurity elastic collisions: 



tt/I 

If— — 

^TT- [J{k',k,x,t) - J{k,k',x,t)]dk', (5) 
27r J 

—rr/l 

J{k, k', X, t) = V, - f{x, k, t), (6) 

provided we use the tight-binding miniband dispersion relation, E{k) = A (1 — 
cos kl)/2 (A is the miniband width), and ignore transversal degrees of freedom. 
For simplicity, z/g and z/j will be fixed constants. The electron velocity in Eq. 
(1) isv{k) = E'{k)/h = Alsmkl/{2h). 

To include shot-noise effects, we assume that the particle fluxes J{k, k', x, t) in 
the clastic collision term Qi[f] are replaced by fluctuating fluxes, J{k, k', x, t) = 
J{k, k', X, t) + dJ{k, k', x, t), where (6 J) = 0. Two such fluxes are independent 
elementary processes unless their arguments are identical. For the same pro- 
cess, the correlations are those of a Poisson process, which yields [8]: 



{SJ{ki, k2, X, t)dJ{k[, k'2, x', t')) 
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in which A is the area of the SL cross section. Inserting the fluctuating terms 
in Eq. (1), we obtain a Bohzmann-Langevin equation (BLE) whose right side 
contains an additional term ^{x, k, t), with zero mean and correlation 



t)e(x , t')) = 6{x - x')5{t ~ t') G{k, k', X, t), (8) 

TTZ/' 

G{k, k', X, t) = — ^ [f{x, k, t) + f{x, -k, t)] [d{k - k') - 5{k + k')]. (9) 



We shall now derive a reduced balance equation for the electric fleld. In order 
to do this, we shall assume that the electric field contribution in the BLE is 
comparable to the deterministic collision terms and that these terms dominate 
the other three. This is the so-called hyperbolic limit, in which the ratio of 
df /dt or v{k)df /dx to (eF/H) df/dk is of order e ^ 1 [4]. Let vm and Fm 
be electron velocity and field scales typical of the macroscopic phenomena 
described by the sought balance equation; for example, let them be the positive 
values at which the (zeroth order) drift velocity reaches its maximum. In the 
hyperbolic limit, the time to it takes an electron with speed vm to traverse a 
distance xq — eFMl/ieNjj), over which the field variation is of order Fm, is 
much longer than the mean free time between collisions, z/~^ ~ h/(eFMl) = ti. 
We therefore define e = ti/to = ^vmNd/{bFmP) and formally multiply the 
fluctuating term ^ and the two flrst terms on the left side of (1) by e. After 
obtaining the number of desired terms, we set e = 1. The solution of Eq. (1) 
for e = is straightforwardly calculated in terms of its Fourier coefficients as 



f^'\k;n)= Y: ffe^''' (10) 

j=-oo 



Fm el y i/e 

Note that Eq. (3) imphes /o^°^ = = n. 

Once /(°) is known as a function of the unknown electron density and electric 
field, we make the Chapman- Enskog ansatz [4] : 



oo 

fix, k, t; e) = /(°) {k; n) + ^ /("^) (k; n) (13) 

m=l 

^^Y.N^'^\n)e-. (14) 

'^^ m=0 

The coefficients f^"^\k; n) depend on the 'slow variables' x and t only through 
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their dependence on the electron density and the electric field (which is itself 
a functional of n). The electron density obeys a reduced evolution equation 
(14) in which the functionals N^"^^(n) are chosen so that the f^"^\k;n) are 
bounded and 27r/Z-periodic in k. Moreover the condition, 



J f^"'\k;n) dk = 2TTft^/l = 0, m > 1, 



ensures that f^"^\ m > 1, do not contain contributions proportional to the 
zero-order term f^^\ N^'^\n) can be found by integrating (1) over /c, using 
(3), and inserting (13) in the result: 



7VM(^)^_/^ / (15) 
Then integration of (2) over x yields a form of Ampere's law: 



^-qT + Y^'"' J v{k)f^"^\k;n)dk = J{t), (16) 
where J{t) is the total current density. 

To find the equations for f^™-)^ we insert Eqs. (13) and (14) in (1), and then 
we equate like powers of e. The result is the following hierarchy of linear 
nonhomogeneous equations: 



+ C^'\x,k,t), (17) 
+ ^^'\x,k,t), (18) 





_ ^ f (0) 

dt^ 



and so on. We have defined Lu{k) = eFh~^du{k) / dk + + i^i/'2)u{k) + 
Uiu{—k)/2, and the subscripts and 1 mean that dn/dt is replaced by N'^'^\n) 
and by N^^\n)^ respectively. In Eq. (17), the correlation of the noise i^^\x^ k, t) 
contains only the zeroth order distribution function. The linear equation Lu — 
S has a bounded 27r/Z-periodic solution provided /^^'^^ Sdk — 0. This solvabil- 
ity condition together with Eqs. (17), (18), etc. also yield the previously found 
jVM of Eq. (15) and the reduced equation (13). 

The solution of Eq. (17) is 
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/'■' = Z ^'^'rlr'-'"'' (19) 

in which S*]^^ ('S'o^'' = 0) is the jth Fourier coefficient of the right hand side 
of Eq. (17). Using Eqs. (10), (11) and (19), we explicitly write two terms 
in Eq. (16), thereby obtaining (after setting e — 1) the following stochastic 
generalized drift- diffusion equation (SGDDE) for the electric field: 



dF f^OFXeNof el dF\ ^ f ^ dF\ d^F 
= -6J{x, t) + A (f, ^] Jit), (20) 



dx ^ 

2e vMFl, [Fl - (1 + 2tI) F^i 
el{ve + y^){Fl + F^f 



A ^ 1 + ^; T ' (21) 
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NnJ h{M)lMkBT)' 

D = ^ A_4n^^\ 

S/iere (F2, + F2) 1^' A/ j' ^^'^ 

^ ibFl - AF^)m2 _ AhvMFUFl - F^){t, + r-^)(nm)^ 

(F|, + 4F2)2m A/(F2, + F2)3 ' ^ > 

^ n{Fi,-2F'){nm2y Shv m jl + T!)[FFM{nm)f 

Fl + AF^ ^ Al{Fl + F^)^ ■ ^ ^ 

Here the electron density is given by the Poisson equation (2), 6 = A/(2A;^T) 
and g' means dg/dn. The fluctuating current density 5J{x,t) in Eq. (20) has 
zero mean and correlation 

/X T( Tf ' +'\\ e^A^Uil5{x-x')5{t-t') ( AhreVMnfh2 \ , . 

The effects of shot noise on the solutions of the deterministic equation are more 
noticeable near bifurcation points. Thus they should affect more evidently the 
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Fig. 1. (a) Self-oscillations of the current density far from bifurcation voltages, (b) 

Zoom of the time interval marked in (a) showing the effects of the shot noise. Data 
correspond to a 100-period GaAs/AlAs SL with well and barrier widths of 5.13 and 
0.87 nm, respectively, and a 3D doping density of 1.4 x 10^"^ cm~'^, as in [13]. The 
units of current density and time are Jq = 1.9x 10^ A/cm^, to = 9.7x lO^^^s, respec- 
tively, and the voltage is 115% of the threshold voltage above which self-sustained 
oscillations appear. 

onset of the self- sustained oscillations of the current that occur in SCSLs [13]. 
These stable time-periodic solutions of the deterministic model are due to pe- 
riodic charge dipole motion from one contact to the other and their frequency 
agrees with experimental observations [4]. Far from the bifurcation voltages, 
the effect of shot noise is small, as shown in Fig. 1. However, near the voltage 
value at which self-oscillations start in the deterministic model, the effects of 
shot noise are more noticeable. Fig. 2 shows that the shot noise may induce 
self-oscillations of the current at voltage values below threshold, for which 
the deterministic model shows relaxation towards the stationary current. A 
great enhancement in the variance of fluctuations near threshold should also 
occur, similarly to the effects of noise on the onset of Gunn oscillations in bulk 
semiconductors [5], cf. Section 7.6 of [9]. 



3 Shot-noise effects in weakly-coupled superlattices 



Deterministic charge transport in weakly coupled SLs is described by discrete 
balance equations [2]. To this day, no one has derived these equations from 
quantum kinetic theory, and therefore we cannot study the effects of shot 
noise on discrete balance equations by perturbative methods similar to those 
outlined in Section 2. What we can do is to add fluctuating particle fluxes to 
the discrete balance equations and assume Poissonian statistics for them. 

Let us consider a WCSL with -|- 1 barriers and wells. The zeroth bar- 
rier separates the injecting region from the first SL well, whereas the A^th 
barrier separates the A^th well from the collecting region. Assume that Fi is 
the average electric field across one SL period, consisting of the ith well and 
the {i — l)th barrier. Similarly, rii is the 2D electron density at the ith well. 
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Fig. 2. Sclf-oscillations of the current density induced by the shot noise for volt- 
ages below the threshold voltage, (a) and (b) show the oscillatory relaxation of 
the current towards its stationary value in the absence of noise, (c) and (d) show 
the self-sustained oscillations induced when the shot noise is present for the same 
voltage value. Numerical values arc as in Fig. 1, but the voltage is 99.3% of the 
threshold voltage above which self-sustained oscillations appear. 

concentrated in a plane perpendicular to the growth direction inside the ith 
well. Then, the Poisson equation (averaged over the ith period) and the charge 
continuity equation are 



Fi-F,_^ = ^{ni-No). (31) 

e = Ji-i^i — Ji^i+i- (32) 

Here Jj^i+i is the tunneling current density across the ith barrier, i. e., from 
well i to well i + 1 with i = 1, . . . ,N. Assuming that the intersubband scatter- 
ing times are much smaller than the escape time from a well, which, in turn, 
is much smaller than the dielectric relaxation time, the main contributions 
to Ji-^i+i are due to sequential resonant tunneling. Electrons in the lowest 
subband Ei of the ith well tunnel to one of the excited levels of the {i + l)th 
well, and then immediately relax to the lowest level thereof. In the limit of 
small (broadening of subbands due to scattering), the following expression 
for the deteministic tunnehng current, Jj^j+i, holds [2]: 



H+l — 
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m*kBT 



In 
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, (33) 
(34) 
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The forward drift velocity v^f\Fi) is a sum of Lorentzians centered at the 
resonant field values Fi = {Ey — Ei)/ {el). Ti{Ei) is proportional to the trans- 
mission coefficient through the ith barrier. The tunneling current is a lin- 
ear function of rii, but it is a strongly nonlinear function of n^+i. Moreover, 
Ji-^i+i ~ ev'^f\Fi)ni/l, for Fi of the order of the first resonant value or larger. 
For such values, the resulting tunneling current density has the same shape as 
assumed in the original discrete drift model [6] . 

Differentiating Eq. (31) with respect to time and inserting the result into 
Eq. (32), we obtain 



dp. 

+ J,^,+i = J(i), (35) 

in which the total current density J{t) is the same function for all i. Equa- 
tion (35) is a discrete version of Ampere's equation. We complete our descrip- 
tion by adding the voltage bias condition 



(AT + i) ^/^-^N + i)r ^''^ 

for the known voltage V{t). 

Our discrete system of deterministic equations consists of Eq. (31) for i = 
1, . . . , iV, Eq. (35) for i = 0, . . . , A^, and Eq. (36). In total, we have (2A^ + 2) 
equations for the unknowns rii, . . . , n^, Fq, . . . , F^, and J{t). We need to 
specify the constitutive relations Jq^i (tunneling from the injecting region to 
the SL), and Jn^n+i (tunneling from the SL to the collecting region). Equa- 
tion (35) evaluated for i = and i = N determines the boundary conditions. 
As tunneling currents at the boundaries, we use the following phenomenolog- 
icai expressions [12]: 



Jo-^i — (^Fq, Jtv-^-w+i — — ^7 ■ (37) 

These conditions are particular cases of the tunneling currents described in 
[2]. The deterministic model consists of Eqs. (31) and (33) - (37). 

In the stochastic model, zero-mean random currents 5Ji_»j+i(t) have to be 
added to the tunneling current densities Ji_>i+i: 

Ji-^i+l — Ji-*i+l + ^Ji^i+l{t)- (38) 

The currents 5Ji_>i+i have correlations [10] 
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for i — 1, . . . , N — 1. The idea behind this form of random tunnehng current is 

that uncorrelated electrons are arriving at the i-th barrier with a distribution 
function of time intervals between arrival times that is Poissonian [8]. More- 
over, the correlation time is of the same order as the tunneling time, so that 
it is negligible on the longer time scale of dielectric relaxation. 

The high-temperature limit of the stochastic model has been numerically 
solved in [10]. Here we shall study a particular effect of shot noise on the 
nonlinear dynamics of a WCSL. The current-voltage characteristics of the cor- 
responding deterministic model consist of a number of stable solution branches 
corresponding to electric field profiles displaying two domains separated by a 
domain wall. The current-voltage plane corresponds to the usual bifurcation 
diagram of a norm of a solution (the total current density) versus the bifur- 
cation parameter (the voltage). We should note that two neighboring stable 
solution branches are typically connected via an intermediate unstable solu- 
tion branch at two limit points that are saddle-node bifurcations. If the initial 
state corresponds to a voltage on the middle of a stable branch, and is 
suddenly switched to a final voltage corresponding to the next branch, Vf, the 
domains relocate after a certain delay time, called relocation time [12]. 

The possible scalings of this mean relocation time can be calculated by means 
of bifurcation theory and the classical results for escape of a Brownian particle 
from a potential well. The idea is as follows. Let (Vl, Jl) be the limit point 
connecting one stable and one unstable branch of deterministic stationary 
solutions via a saddle-node bifurcation: two solutions exist for V < Vl, and 
they disappear for V > Vl- Let us assume that Vi < Vf < Vl, and that for 
V = Vf three stationary solutions coexist with current densities Ji < J2 < J3. 
The solution branch corresponding J = J3 is the same stable branch containing 
the initial point. This branch coalesces with the branch of unstable solutions 
corresponding to J2 at the limit point (Vl, Jl)- Provided (Vf — Vl) is of the 
same order as the noise amplitude in the appropriate dimensionless units, a 
projection of the noise on the eigenvector corresponding to the zero eigenvalue 
in the linear stability problem at Vf = Vl appears in the amplitude equation 
for the saddle-node bifurcation. The latter becomes the following stochastic 
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amplitude equation: 



da 
'di 



a{Vf-VL)+Pa' + ^m. 



(42) 



In this expression, a is proportional to \J—Jl\ and ^(t) is the zero-mean, delta- 
correlated white noise. The positive parameters a, j3, and 7 can be calculated 
by projection methods [14], but their precise form does not matter for the 
argument we want to make. Note that Eq. (42) has two stationary solutions 
with \ J-Jl\<x \Vf - Vl|i/2 if 7 = 0. 

For Vf < Vl, Eq. (42) describes the escape of a Brownian particle from a 
potential well corresponding to the cubic potential C/(a; V/ — Vl) = ol (V/ — 
V/,) a + j3a^ /?). Provided the height of the barrier is large compared to the 
noise strength, the reciprocal of the mean escape time is proportional to the 
equilibrium probability density P — e~^^'^ /Z evaluated at the maximum of 
the potential. This yields 



cf. Eq. (8.6.17) in [15]. Thus, there exists a relatively large voltage interval 

a 3/2, over which the logarithm of the relocation time 
scales superlinearly with | V}— Vl|. In terms of the current, ln(Tre/oc) ^ \ J—Jl\^-i 
because a(x\J ~Jl\ = 0(|V} - Vl\^I^). 

There is some evidence that the parabolic region near the limit point might be 
so narrow that the amplitude equation (42) breaks down. For example, Fig. 5 
of Ref. [16] shows that the maximum and the limit point of the current-voltage 
characteristic are extremely close for all branches. In this case, it may very 
well happen that the region -f'^/^ (3^^^ / a < \Vf -Vl\ <^ K in which Eq. (42) 
holds is too narrow. This case might be better approximated by analyzing the 
effect of shot noise on a spiky limit point, at which the slope of the current- 
voltage characteristic is discontinuous, and the term proportional to in (42) 
is replaced by a piecewise linear function of a. Then a simple calculation of the 
barrier height yields \n{Treio<^ oc |V/ — Vi,p oc \ J — Jl\^ instead of Eq. (43). The 
numerical solution of the stochastic equations could be used to discriminate 
between this prediction and Eq. (43). There are early experiments in which 
\n{Treioc) was fitted to a straight line [11], but no theoretical prediction existed 
at that time and perhaps more careful measurements should be attempted in 
order to validate the theory. 




(43) 
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4 Conclusions 

We have modelled the effects of shot noise in strongly and weakly coupled 
superlattices. For a SCSL, we have proposed a Boltzmann-Langevin equation 
containing fluctuating terms that represent shot noise. Then we have used 
a consistent Chapman-Enskog method to derive a stochastic drift-diffusion 
equation for the electric field. In the case of WCSLs, we have slightly general- 
ized the discrete Langevin equations proposed in [10], and found the scaling 
for the mean relocation time that electric field domains take to react to a 
switch in the applied voltage. 



The figures in this paper were calculated by Guido Dell'Acqua and Ramon 
Escobedo, to whom 1 am very much indebted. This work has been supported 
by the MCyT grant BFM2002-04127-C02, and by the European Union under 
grant HPRN-CT-2002-00282. 
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